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Abstract 

We consider the non-commutative inflation model of [3[ in which it is the unconventional disper- 
sion relation for regular radiation which drives the accelerated expansion of space. In this model, 
we study the evolution of linear cosmological perturbations through the transition between the 
phase of accelerated expansion and the regular radiation-dominated phase of Standard Cosmology, 
the transition which is analogous to the reheating period in scalar field-driven models of inflation. If 
matter consists of only a single non-commutative radiation fluid, then the curvature perturbations 
are constant on super-Hubble scales. On the other hand, if we include additional matter fields 
which oscillate during the transition period, e.g. scalar moduli fields, then there can be parametric 
amplification of the amplitude of the curvature perturbations. We demonstrate this explicitly by 
numerically solving the full system of perturbation equations in the case where matter consists of 
both the non-commutative radiation field and a light scalar field which undergoes oscillations. Our 
model is an example where the parametric resonance of the curvature fluctuations is driven by the 
oscillations not of the inflaton field, but of the entropy mode. 
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I. INTRODUCTION 



Non-commutativity of space and time is a general prediction for an effective field theory 
coming from string theory pj . This leads to a maximum wavenumber pc for any field which 
is subject to this uncertainty. In turn, this leads to a dispersion relation which is modified 
at high clensities. In particular, regular radiation is described by such a modified dispersion 



relation [2]. As shown recently in 3|], the dispersion relation will likely contain two branches 
(two frequencies for any wavenumber smaller than the maximal momentum). This two- 
branch structure may be related to the T-duality symmetry of string theory which has been 
made use of in a different context in string cosmology in {4]. At high temperatures, both 
branches of the dispersion relation are occupied, whereas at low temperatures only the lower 
branch, the branch which for p <^ Pc has the regular linear form, is occupied. 

In [3| it was realized that for suitable classes of dispersion relations, regular radiation 
leads to inflationary expansion of space at high temperatures. The reason is that as space 
expands and the physical momentum of a flxed comoving mode decreases, the energy of 
the mode increases. This increasing energy leads (in the context of a background described 
by Einstein gravity) to accelerated expansion. We called this model "non-commutative 
inflation" . 

Recently we have studied the generation of cosmological fluctuations in non- 
commutative inflation. In contrast to simple scalar fleld-driven inflationary models, the 
fluctuations are thermal of origin and not quantum vacuum fluctuations, the reason being 
that matter in our scenario is a thermal bath rather than an ultracold Bose condensate of a 
scalar fleld. In the limit that the accelerated expansion is almost exponential, the spectrum 
of fluctuations is almost scale-invariant, as should be expected from the general symmetry 
arguments of 6|]. 

An interesting aspect of non-commutative inflation is that the transition between the 
period of accelerated expansion and the radiation phase of standard cosmology is straight- 
forward: matter always is in regular radiation. As the universe expands, the equation of 
state of this radiation changes from being that of an accelerated universe early on to that of 
ordinary radiation. Thus, there is no need for any period of "reheating" , a rather non-trivial 
phase in usual scalar fleld-driven inflationary models. 

In this paper we study the evolution of the fluctuations through the transition between 



the inflationary phase and the phase when the equation of state changes to that of ordinary 
radiation. Since there are no oscillations in the equation of state of that matter giving rise 
to inflation (in scalar field-driven inflation the phase of accelerated expansion ends with 
oscillations of the scalar field giving rise to inflation which translates into oscillations in the 
equation of state of matter). As a consequence, there is no parametric resonance instability 
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of matter [7, 8 or metric 
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121 ] fluctuations in this phase, at least in the absence of 



other matter fields which undergo oscillations. 

In models of matter beyond the Standard Model, in particular in models motivated by 
superstring theory, there are many scalar matter fields which are expected to undergo oscil- 
lations at the end or after the period of inflation. These fields lead to entropy modes. Here, 
we will study the evolution of cosmological fluctuations in the presence of such oscillating 
entropy modes (oscillating "moduli" fields). We find that these entropy modes can indeed 
lead to parametric resonance of the curvature fluctuations on super-Hubble scales. 

The importance of parametric resonance instabilities in the process of reheating at the 
end of inflation was first recognized in j3] (see also Q]). The resonance of matter fluctuations 
during reheating in scalar field-driven inflation was then studied in more detail in [isl (where 
the name "preheating" for this resonance phenomenon was coined) and [l^. A detailed 
discussion of the efficient broad-band and the much less efficient narrow-band resonances 
was given in ISj. It was ffist conjectured in 9|] that super-Hubble (but sub-horizon) scales 
metric fluctuations may also be parametrically amplified. In particular, there are no causality 
constraints |l^ which prohibit this (recall that in inflationary cosmology the horizon is larger 
than the Hubble radius by a factor of exp(7Vp), where N^. is the number of e-foldings of the 
inflationary phase). However, as shown in [101], in the case of purely adiabatic fluctuations 
there is no resonance at linear order in perturbation theory. However, it was then established 
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121 ] that oscillations in the equation of state of matter can lead to an increase of the 



curvature fluctuations on super-Hubble scales in the presence of low mass (mass smaller 
than the Hubble rate during inflation) entropy modes. 

In previous work, most of the attention has focused on oscillations of the inflaton driv- 
ing the parametric resonance of the entropy modes of cosmological fluctuations, which in 
turn feeds the exponential growth of the curvature perturbation via the usual sourcing of 
curvature fluctuations via an entropy inhomogeneity. In the context of non- commutative 
inflation, it is the background value of the entropy field which undergoes oscillations. Here, 



we show that such oscillations may also be able to excite a parametric resonance for entropy 
fluctuations. 

The outline of this article is as follows: In the following section we give a brief review 
of non-commutative inflation. Then, in Section 3, we study the growth of cosmological 
fluctuations during the transition between accelerated expansion and the usual radiation 
phase in the case of a single component of matter, the radiation. As expected, we find 
constant curvature fluctuation on super-Hubble scales. The key section of our paper is 
Section 4, in which we study the equations of motion for cosmological fluctuations for non- 
commutative radiation coupled to an oscillating scalar field (the modulus field acting as the 
entropy mode). We derive the relevant perturbation equations, solve them numerically, and 
give some analytic insight into the solutions. We conclude with a discussion of our results. 



II. NONCOMMUTATIVE INFLATION 

The starting point of the non- commutative inflation model of [3| is the modified dispersion 
relation for massless particles 

-p^c^fiEf = (1) 
which results from the non-commutativity of space and time. Here, 

fiE) = l + {\Er. (2) 

Note that p and E denote momentum and energy, respectively, a > 1 is a positive constant, 
and the length scale A determines the maximal momentum pc- For a > 1 the dispersion 
relation has two branches. On the upper branch, the energy increases as the momentum 
decreases while the universe is expanding. It is this behavior which, for a range of values 

n 

of a explored in [3|], leads to inflationary expansion of space when the usual Friedmann- 
Robert son- Walker equations for the coupling of the background space-time to matter are 
used. 

The modification of the dispersion relation leads to a deformed thermal spectrum 
The energy density p is given by 
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fE 
f 



(3) 



and the expression for the pressure V is 

lfpiE)dE 

In the high energy hmit, XE ^ 1, f'E/f ^ a and p (x T. This approximation leads to a 
nearly constant equation of state parameter w = V/p. 

Regular commutative inflation has no intrinsic entropy fluctuations since 

where 6 indicates spatial variations. The same result holds for non-commutative radiation. 
This can be seen by taking on one hand the time derivatives of Eqs. ([3]) and (jl]), and on 
the other hand the infinitesimal spatial variations of these quantities. The vanishing of the 
intrinsic entropy fluctuations will be important later on. 

If the background geometry is described by the Einstein action, the Friedmann and energy 
conservation equations take their usual form 



p = -3H{l + w)p. (7) 

Inserting the expressions and (jl]) into the above equations, it follows that at high temper- 
atures, a period of accelerated expansion results provided that the parameter a is suitably 
chosen ^ (values of a slightly larger than 1). The resulting evolution of the equation of 
state parameter w and of the adiabatic sound speed is shown in Fig. [TJ 

Inflation has to last more than 60 number of e-foldings (this number assumes that the 
energy scale at which inflation takes place is close to its upper bound) to solve the problems 
of Standard Cosmology which inflation was designed to address. If noncommutative inflation 
ends roughly when XT^nd — 1? then the number of e-foldings is given by 

= ln[a{tend)/a{U)] = f"" Hdt = ^ JfML^dT . (8) 

In the high energy limit, we obtain 

1 , XTi InATj 

3(1 + XTend 3(1 +w) 

where we have used XTend ~ 1- 
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FIG. 1: Evolution of the equation of state parameter w and of the adiabatic sound speed 
for a noncommutative fluid during and after the period of inflation. The value of a is 

set to 1.099. 

III. EVOLUTION OF NONCOMMUTATIVE FLUID PERTURBATIONS 
A. Equations of motion 



In this section we will take matter to consist only of the non-commutative fluid. In 
this case, since there are no entropy fluctuations we expect the curvature fluctuations to 
be conserved on super-Hubble scales. In addition, since there is no oscillating matter field, 
there cannot be any parametric resonance effects. 



We will work in longitudinal 



gauge to study the linearized equations of motion for cosmo- 



logical perturbations (see 16|, ll7| for comprehensive reviews of the theory of cosmological 
perturbations and [l^ for a pedagogical overview). In this gauge, the metric takes the form 



ds'' 



2^)dt^ + 0^(1 - 2<iJ)-fijdx'dx^ 



(10) 



where is the spatial part of the background metric which we in the following will take 
to be Euclidean. The perturbations are described by $ (the Bardeen 1^ potential) and 
both functions of space and time. 

If matter consists only of the noncommutative (NC) fluid, the perturbed energy- 



momentum tensor can be written as 

5T\ = -5p, 6T\ = 6q„ 5T, = 5V5], (11) 

where we have assumed vanishing anisotropic stress. 
The perturbed Einstein equations are 

2,H{H^ + ^) - ^V^^ = -ATiG5p, (12) 
H^ + if = -47rG5g, (13) 
^ + mi) + H^ + {2H + 3iJ2)$ _ ^V2(* - $) = AtxG5V. (14) 

The vanishing of the anisotropic stress tensor makes it possible to set $ = \[^. 

The conservation of energy-momentum tensor for the NC fluid leads to the following 
equations for 5p and 5q, 

6p + 3H{6p + 6V) - 3p(l + w)^ + \v^Sq = 0, (15) 
5q + 3HSq + p(l + w)^ + Sp = 0. (16) 

Combining Eqs. ( fT2l) and iHM and working in momentum space ^ —k'^^k, we can 
obtain a single differential equation 

+ (4 + 3c2)i7$fc + + + 3/7^(1 + $fc = , (17) 

where we have used 6V = c^Sp. 

For adiabatic fluctuations, and in particular in our single perfect fluid system, = c^. 
However, in non-adiabatic fluids, for example systems with several fluids or scalar fields, 
there may exist intrinsic isocurvature perturbations which are denoted by 

pV = (cl - cl)6p. (18) 



The Sasaki-Mukhanov variable 



20 



21| for a fluid in terms of which the action for cos- 



mological perturbations has canonical form, and which is thus useful for the quantization of 



these fluctuations, can be written as f = ac^ ^Qf where 17|, |22| 



Qf = ^^i^^-^-^"^]- (19) 



Vp + v V H 

If matter is a scalar field, this gauge-invariant quantity becomes 



Q^ = 5<P-j^^. (20) 



In usual scalar field-driven inflationary models, the scalar field begins to oscillate about 
the minimum of its potential after the end of inflation. These oscillations lead to a singularity 
if Eq. (fT7|) is used because of the 1/0 terms in the coefficients. To avoid this singul arity 



during the preheating phase, the Sasaki- Mukhanov variables must be used instead of $ 



23|. 



Even in the absence of oscillations, if the change in the equation of state of the background 



in smooth, it is safer to use the Sasaki-Mukhanov variables instead of $ (see e.g. 2J, 1251]). 

The situation, however, is different in our fluid-driven inflationary model. The speed of 
sounds Cg is negative during the inflationary period. Then, after the end of inflation, 
crosses zero and evolves to its final positive value cj. = 1/3 (see Fig. [1]). The evolution 
equation for Qf, Eq. (\A2\i . contains (c^)~^ terms which diverge when = 0. Therefore it is 
convenient to use the $ equation of motion ( fTTl) to evolve the fluctuations in our fluid-driven 
inflation model in order to avoid the singularity at = 0. 

The curvature perturbation ( on uniform energy density hyper-surfaces is related to the 
Sasaki-Mukhanov variable as 

( = $ - :^($ + = — T^=Qf. (21) 

The evolution equation for is 

a = f^$.-^T^pr. (22) 

H a p(l + w) 

In the case of adiabatic fluids when [pT = 0), then for long- wavelength perturbations (k is 
conserved. 



B. Numerical results 

We performed numerical calculations of the evolution of the linear cosmological pertur- 
bations during and after NC inflation. The simulations were done using the value k = 5aiHi 
for the comoving momentum (the subscripts i standing for the initial time of the analysis) 
and started from ATj = 9900 to satisfy the condition of Eq. (I9l). Thus, the modes are inside 
the Hubble radius in the initial stage of the inflationary period. Quantum vacuum initial 



conditions for the fluctuations were used 



30l ]. To implement these initial conditions, we 



postulate the usual harmonic oscillator vacuum initial conditions for the Fourier modes of 
the variable v in terms of which the action for fluctuations has canonical kinetic term. This 
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FIG. 2: Evolution of the Sasaki-Mukhanov variable Qfk in the case when matter is a single 
NC fluid. Also shown are the curvature perturbation and the Bardeen potential 
The simulations starts during the inflationary period and the results shown are for the 

values = 10~^mpi and a = 1.099. 

variable is given by 



where the variable z is 



Vk = —Qfk = -zCk, 



The variable Vk obeys the following wave equation 



(23) 



(24) 



v'k+{k'cl--]vk = 0, 



where a prime denotes the derivative with respect to conformal time t] {dt 
initial conditions for the A;'th Fourier mode of v are given by 



(25) 
adrj). The 



Then, from Eqs. fl^T]) and fl22]) . we obtain the initial conditions for $ (suppressing the index 
k for now) 



where Vi = v{ti). 

We set the cutoff mass A^^ = IQ'^nipi to be consistent with the observational data 
and use a = 1.099, a value which gives an inflationary background. With the initial 
conditions described above and these parameter values, we numerically solve the Eq. ( ITTj) . 

As shown in Section 2, the speed of sound is given by = c^. Thus, during the inflation- 
ary phase is negative. This leads to an instability of short wavelength fluctuations. Since 
the focus of this paper is on the evolution of the fluctuations on super-Hubble scales, and 
since - according to the prescription of - we should impose initial conditions for the 

fluctuations at the typical thermal wavelength, and modes thus do not spend a long time 
on sub-Hubble scales, we will not worry about this instability in this paper. 

In Fig. [21 we display the evolution of ^k.Qfk and (k during and after inflation. Fig. [2] 
shows that the curvature perturbation stays constant for super-Hubble scales both during 
inflation and after the end of inflation. This is the same behavior as holds for adiabatic fluc- 
tuations in scalar field-driven inflationary models. Meanwhile, the time-dependent equation 
of state causes the Bardeen potential and also Sasaki-Mukhanov variable Q fk to increase 
for super-Hubble scales during inflation as time increases. We wish to emphasize that, as 
expected, there is no amplification of the amplitude of curvature perturbations at the end 
of infiation. This can be understood from the fact that in Eq. (!25|) . the function z" /z term 
does not include any oscillating terms at the end of infiation. 

Our numerical analysis (Fig. [2]) also shows that the curvature perturbation (k as well as 
<l>fc and Qfk increase on sub-Hubble scales . This is unlike what happens in the usual scalar 
field-driven infiation model. The difference, as discussed above, is due to the fact that in 
our case the sound speed is negative during inflation. 




(27) 




(28) 



IV. EVOLUTION OF NONCOMMUTATIVE FLUID PERTURBATIONS IN- 
CLUDING AN OSCILLATING SCALAR FIELD 



Particle physics models beyond the Standard Model, in particular those based on super- 
symmetry, supergravity and superstring theory, typically contain many scalar fields. These 
fields will be displaced from the minima of their potentials by quantum fiuctuations during 
the period of infiation. After infiation, they will start to oscillate about their minima. In 
this section we will study the effects of these oscillations on the parametric excitation of 
entropy fluctuations on super-Hubble scales. 

In different contexts, cosmological perturbations in models with both scalar flelds and 
fluids have been studied before. For example, in warm inflation models 



the scalar 



inflaton fleld is coupled to a thermal fluid. Since the scalar fleld evolution is typically 
overdamped, there will be no parametric resonance effects. The growth of fluctuations in 
this context has been considered e.g. in 27j. Our situation is different in that the dynamics 
is driven by the fluid rather than by the scalar fleld. In a more general context, perturbations 
in a system consisting of a scalar fleld plus an ideal gas were recently also considered in 28| , 
in particular with applications to quintessence cosmology in mind. 

In the following, we flrst discuss the equations of motion for the cosmological fluctuations, 
and then present our numerical results. 



A. Equations of motion 

The background scalar fleld cf) satisfles the Klein-Gordon equation 

+ 3i^0 + = 0, (29) 

where = dV/dcl). In the presence of both a fluid and a scalar fleld, the background 
Friedmann equation is 

H^=^-^ [pf + + V{<P)^ . (30) 

where subscript / denotes a fluid quantity. 

In Fig. [3l we plot the evolution of a free scalar fleld in a universe dominated by 
non-commutative radiation. Different plots are for different values of the scalar fleld mass 
m. 
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FIG. 3: Evolution of the scalar field during and after the period of NC inflation for various 



values of the scalar field mass m. The parameters A ^ = 10 ^m„i and a = 1.099 were used. 



The total perturbed energy density, pressure and momentum are 

6p = 6pf + 6p^, 6V = 6Vf + 6Vf, 
6q = 6qf - 



where 



(31) 

(32) 
(33) 



By inserting the metric including cosmological fluctuations into the curved space-time 
Klein-Gordon equations, one obtains the following equation for the fluctuation of the scalar 
field 



4$i 



(34) 



By subtracting times the 0—0 perturbed Einstein equation from the space-space perturbed 
Einstein equation, one obtains the following evolution equation for the Bardeen potential 



+ (4 + 3ci)H<^, 



2H + 3H'{1 + ci) + 47rG(l - c; 
47rG[(l - cl)<p6<Pk - (1 + clWM, 



2\ J,2 



(35) 



where 



H=-A7iGipf{l + Wf)+<i)^). (36) 

Note that in the above is the speed of sound of the fluid component. 

The isocurvature perturbation for a multi-component system (the components being la- 



belled by indices ...) takes the form 16| 



where we have used the abbreviations 

h = p + V, hi = p^ + Vi. (38) 

The first term in Eq. ( 1371) represents the non-adiabatic pressures of the individual com- 
ponents and the second term the relative isocurvature perturbations between the different 
components. If energy transfer between the components is neglected, Sij becomes 

= - (39) 

In our case we consider one fluid and one scalar field, and the relevant term, the expression 
for Sfcj, is 

^ 3^^$^ + 3if$,. + g'I'fc pf{l + wf) + 0^ 
47rGp;(l+^,^) Pf{l + Wf)^^ 



AnG 1 - c 



cl^/Pfi L+W] 



2 



H cl 

Then, from Eq. (1371) . the relative isocurvature perturbation is 



Pf{l + w)Qf). (40) 



2 pf{l + Wf) + 02 
where the speed of sound of the scalar field component is 



c^, = ^ = l + ^- (42) 
^ P4> 3H(j) ^ 



As we have discussed in Section 2, there is no non-adiabatic pressure for the NC fluid. The 
non-adiabatic pressure for the scalar field can be written as 

P^int = (43) 
2Vs 



+ w<^) + - v^p/(i + ti^/)g/). 



2,H(i) P 
The sum of Eqs. (1411) and (H3|) gives the total non-adiabatic pressure which is the source 
term in the equation of motion for the curvature perturbation ( (see Eq. fl22l) ). 

In the case of a mult i- component system, the variable ( which geometrically represents 
the curvature perturbation on the uniform energy density slices is given by ^| 

^ -\/Pf{'^ + Wf)Qf + ^Q4> 
and its equation of motion is given by the generalization of (l22l) as 



= 5^^$,, ^pV (45) 

From this, we see that in the case of the NC fluid coupled to a scalar field, there is a non- 
vanishing source for Ck which is not suppressed on large scales (the second term on the right 
hand side of Eq. (HS]) ). and if the scalar field is oscillating, there is a possibility of parametric 
resonance enhancement of Cfc- We will now turn to a numerical study of this question. 



B. Numerical Calculations 

We have evolved Eqs. and (l35l) numerically for a massive scalar field potential 

V{ct>) = \m''<p\ (46) 

The simulations began when the perturbation modes are still inside the horizon, and the 
scalar field fluctuations were started in their vacuum state given by 

5(l>u{ti) = (50,(t,) = (47) 

Since we are considering initial conditions for which the energy density is dominated by the 
NC fluid, i.e. pf ^ during the early stages of the inflation, we used the same initial 
conditions for $ as in Eq. (l28l) . 
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FIG. 4: Evolution of the individual Q variables, namely Qfk for the NC fluid and Q^k for 
the scalar field, the Bardeen potential and curvature perturbation Qki in our model of 
NC inflation. The results are for the choice of the NC inflation parameters A"^ = 10~^mpi 
and a = 1.099. In the scalar field sector, we set the initial value of (p as (pi = O.lnipi, 

<j)i = 0, and mVif2 ^ iq-2_ 






In Fig. [3] we plot the evolution of the background scalar field during and after infiation for 
different values of the mass parameter of the scalar field. We used the NC fiuid parameters 
A"^ = lO'^nipi and a = 1.099, and chose the initial value of the background scalar field as 
(pi = lO'^rUpi with vanishing initial velocity. We considered different values of the mass. 

For m > H, the scalar field oscillations are rapidly damped during infiation, and no 
oscillations arise during the post-infiationary era. For m <^ H, the scalar field is over- 
damped during and immediately after infiation and thus hardly rolls at all during the time 
interval of our simulations. The interesting case for our purpose arises when the mass m 
is slightly but not much smaller than H. In this case, the scalar field rolls slowly during 
infiation and, when the infiationary phase ends and the Hubble parameter starts to decrease, 
the slowly decreasing Hubble parameter becomes comparable in value to the mass of the 
scalar field and then begins to oscillate about the minimum of its potential. The time scale 
of the oscillations becomes smaller than the Hubble time as H decreases. This provides the 



12|. 



conditions where parametric resonance excitation of the entropy modes is expected 

The evolution of Qf,Q^,^k and (k are plotted in Fig. |4] for = lO'^H"^ and 0j = 
10~^mpi and for the same parameter values of A anda as in Fig. [3l With these parameter 
values, parametric amplification of the curvature perturbation results after the end of the 



inflationary phase. 

To see the region in which parametric resonance occurs in detail, in Fig. \5\ we focus 
on the time interval where the scalar field oscillations occur. For a value of the mass 
given by = 10~^iJ^, the oscillation of the background scalar field during this time 
interval are seen (see Fig. [3]) to cause oscillations of the scalar field fluctuation variable Q^. 
Via the gravitational coupling of the two matter components, this leads to the parametric 
amplification of the fluid fluctuation variable Qf. 

To understand the parametric resonance after inflation in our NC fluid model of inflation 
with an oscillating scalar matter field, we must consider the coupled system of evolution 
equations for Qj and Qif, which are given in the Appendix (Eqs. (jAl]) and ([M])). We 
assume that the time scale of the resonance is short compared to the Hubble time H~^. In 
this case, we can neglect that Hubble damping terms in the evolution equations. We will 
also drop the other terms involving only first derivatives of the Q variables - we later check 
the self-consistency of this approximation - and we will also drop coefficients of Qf and 
terms which are obviously suppressed compared to other coefficients. Then, we obtain 
approximate equations 



Qfk + 

where 



+ GnGApf - 6nGB^' 



Qfk-A7TGJpf{l + Wf)[C + 3(t>D]Q^k ^ 0, (48) 



and 



where 



+ 



A 
B 

G 
D 



[1 + w)il + 3c^) + (1 - c^)(l + c^) - 2(1 + WffPl., 



;i-a + (i+^)(i+a^, 

p 



(1 
(1- 



1^, 



P 



Q^k + SttGJ pf{l + Wf) 



H ^ 



3 + 



3 (1 - c's? _ Hl+cl){l + Wf) pf 

4 c2 4 c2 p' 



(49) 
(50) 

(51) 
(52) 



Qfk^O, (53) 



(54) 



Note that at the end of the inflationary phase, Pf/p ^ 1/2 because of pf ~ p^. This system 
of equations represents two coupled harmonic oscillators with mass terms which, during the 



time interval when (f) is oscillating, contain some periodically varying coefficients. We can 
diagonalize this system of equations to get two decoupled equations of Mathieu type an 



[A{k) - 2qcos2z\Qk = 0, (55) 



where Q stands for one of the eigenvectors of the diagonalization process and z = mt. 

We need to estimate the values of A{k) and q once begins to oscillate and can be 
approximated as = 0e sin mt where (pe is the amplitude of the scalar field at the end of 
inflation. While for g <^ 1, the resonance occurs only in narrow bands for q > 1, the 
resonance occurs for a broad range of values of k [l^. Our numerical simulations indicate 
that we are in the broad resonance region. Let us for a moment neglect the coupling terms 
between the Qf and Qfp equations. In this case, the value of q in the Qf equation would be 
of the order 

g ~ Gtt f — 1 (56) 
\mpij 

and would be less than 1. However, if we estimate the amplitude q of the periodically varying 
coefficient in the cross coupling term in Eq. (148|) . we find 

q . (57) 

mmpi 

Inserting the value of H at the beginning of inflation, we obtain a number which is of 
order unity. Looking at the equation we find a very similar result: the periodically 
varying coefficient multiplying Q f leads to a value of q given by Eq. (l57j) and hence appears 
to contribute to the sourcing of the broad resonance. The problem with this argument, 
however, is that the value of H at the end of inflation has decreased to a value smaller than 
m, and hence the above estimate of q in Eq. (1571) gives a value smaller than 1. Returning, 
however, to the general equations in the Appendix, we see that there are coefficients which 
contribute to q which are enhanced by and hence, during the period when the equation 
of state transits from that of inflation to that of regular radiation, leads to an enhancement 
of q. 

To summarize the discussion in the above paragraph: We have presented analytical es- 
timates which support the numerical results which imply that the parametric amplification 
of the coupled Qf and system is in the broad resonance domain, as shown in Fig. [51 



Note that the Qf and Q^j, equations have a singularity when cj. = 0. However, we are nu- 
merically solving not the Q equations, but the $ equation, and our evolution equations are 
non-singular. 

The amplification of and Qf also gives rise to a parametric amplification of as can 
be seen in Fig. |5l From Eq. (145|) . it can be seen that the large growth of ( is generated 
through the isocurvature source term. We have compared the magnitudes of the curvature 
perturbation and the isocurvature source term in Fig. O This confirms numerically that 
the amplification of ( occurs once the isocurvature source has become sufficiently large. 

As a further consistency check, we can perform a rough estimate of the duration of the 
period of oscillations until the amplitude of the curvature perturbation can start to grow. 
As soon as the scalar field starts to oscillate, the isocurvature perturbation in Eq. psl) grows 
exponentially as exp(/im(t — tg)) where tg is the time when the oscillation begins and the 
Floquet exponent describes the rate of exponential growth of the instability. For q ^ 1, 
unstable modes grows extremely rapidly with /i ~ 0.2 jl3|. Then after integrating Eq. fHSl) . 



the curvature perturbation becomes 



C(t) C(t) H{te)pT{te> /ggN 

^^m{pf{U){l + Wf{U)) + m) 
The time t when the increase in the curvature perturbation becomes visible compared to 
its initial value is given by = C,{t) — C{te) ~ (9(l)C(^e) at the time when the amplitude 
of the curvature perturbation begins to grow. Then we can estimate the duration of the 
oscillatory period until the curvature perturbation starts to increase substantially compared 
to their initial value 

where we have used pjite) — P4,{te),(p(te) ~ m(j)(te), and Wfitf.) ~ m ~ lOif(te), 

0(te) ~ 10~^mp; and Cite) /pT{te) ~ 10^ at the beginning of oscillation for the parameter 
values in Fig. El With these parameters, Eq. ( l59l) gives 

t-te - — ^ 6-^H-\Q ~ 3H~\Q . (60) 
fim jjL m 

From the numerical calculations (Fig. [H]), (t — te) numeric — A ln(a/ao)-ff~^(te) ~ 3iJ^^(te)- 
This is consistent with the above analytical estimate. 



42 44 46 48 

In a/a^ 

FIG. 7: Evolution of Qfk for the NC fluid, Q^k for the scalar field, the Bardeen potential 
the curvature perturbation (k, and pT during the inflationary period for the parameter 
values = 10~^mp; and a = 1.099. We set 0i = 10~^mp;, 4>i = 0, and m? / H"^ = 10~^. 

We end our discussion by showing the results of a simulation with different values of the 
parameters. In Fig. [TJ we plot Qfk, Q(t>ki Cfc and pV for 0j = 10~^mp; and rn? / H"^ = 10~^. 
Once again we find parametric amplification of the curvature perturbation. 



V. DISCUSSION 



In this paper we have studied the evolution of cosmological fluctuations through the 
transition between inflationary phase and radiation phase in our non-commutative inflation 
model 0, 5]. In this model, there is no traditional reheating phase with oscillating inflaton 
field. Thus, if there is no matter entropy mode which oscillates, there is no possibility for 
a resonance growth of curvature fluctuations on super-Hubble scales. As we show here, if 
a scalar matter field happens to be oscillating at the end of inflation, which is expected to 
happen in many models beyond the Standard Model of particle physics with light moduli 
fields (with mass of the order of the Hubble constant during the inflationary phase), then 
there can be parametric resonance of curvature fluctuations at the end of inflation on super- 
Hubble scales. This resonant growth happens provided that the scalar field mass at the end 
of the period of inflation is comparable to the Hubble rate at that time. 



The resonance we study here is different from the usual resonance in scalar field-driven 
models, in which the resonance is driven by the coherent oscillations of the dominant compo- 
nent of matter. In our situation, it is oscillations of a sub-dominant mode which is sourcing 
the resonance of the curvature fluctuations. We show that there is an entropy mode which 
undergoes resonance, and once the entropy mode has acquired a sufficiently large amplitude, 
a resonant growth of the curvature fluctuations is induced. 

A characteristic feature of our non- commutative inflation model is the fact that the fluid 
speed of sound cj. is negative during the phase of accelerated expansion. This appears to lead 
to an instability of fluctuations on small scales. Since this work focuses on the transition 
of fluctuations through the period when the equation of state of the background changes, 
rather than on the question of initial conditions for the inflationary phase, we here to not 
address this problem. However, we will need to come back to this issue in future work. 
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APPENDIX A: EQUATIONS FOR Qf AND FOR A MODEL CONSISTING 
OF ONE MATTER FLUID AND ONE SCALAR FIELD 



The Sasaki- Mukhanov variables for the fluid and the scalar field are defined in Eqs. (fT9|) 
and fl20|) . respectively. They satisfy the following set of coupled equations: 



Q^k + 3HQ^k - 47rG-^(l - c^^)a/p/(1 + Wf)Qfk + 



16 

+ 247rG'02 



(Al) 



Qfk 



Qfk + 47rG-^(l - 4) + Wf)Q4,k + 



2„2 



k c 



Q7rGpf{l + Wf){l + 3clf) 
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- IQ^'G'P^Pfil + Wf){l + cl,) - ^-Hclf - ^-Hil - + AnG 



p(l + Wf) 



H 



-sf^ 



Qfk 



-AnG^Pfil + Wf) 

1 a2 -I 



[l + C 



sf) 



H 



^vrG-^p^l + Wf) + 30(1 + c%) - 47rG-|^(l + c^^) 



+ 



H6\ 



Q(f>k — 0, 



(A2) 



where we have used c^j = c^j. Similar equations were derived in [27| for constant Wf 



-af - "-sf 



sf - "-af 

c^r and with a friction term. Note that Eq. flA2l) becomes a single fluid equation for 



0. And Eq. (lAip agrees with the equation for a single scalar field equation when pf = 0. 
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